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Abstract 

We analyze the low energy spectrum of bound states of the N=l SU(2) SUSY 
Yang-Mills Theory (SYM). This work continues the investigation of the non-pertur- 
bative properties of SYM by Monte Carlo simulations in the Wilson discretization 
with dynamical gluinos. The dynamics of the gluinos is included by the Two-Step 
Multi-Bosonic Algorithm (TSMB) for dynamical fermions. A new set of configura- 
tions has been generated on a 16 3 • 32 lattice at (5 = 2.3 and k = 0.194. The analysis 
also includes sets of configurations previously generated on a smaller (12 3 • 24) 
lattice at ac = 0.1925,0.194 and 0.1955. Guided by predictions from low energy 
Lagrangians, we consider spin-1/2, scalar and pseudoscalar particles. The spectrum 
of SYM is a challenging subject of investigation because of the extremely noisy cor- 
relators. In particular, meson-like correlators contain disconnected contributions. 
The larger time-extention of the 16 3 • 32 lattice allows to observe two-state signals in 
the effective mass. Finite- volume effects are monitored by comparing results from 
the two lattice sizes. 
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1 Introduction 



The N=l SU(N C ) SUSY Yang-Mills (SYM) theory is the simplest instance of a SUSY 
gauge theory and presently the only one viable for large scale numerical investigations. It 
describes N 2 -l gluons accompanied by an equal number of fermionic partners (gluinos) in 
the same (adjoint) representation of the color group. Veneziano and Yankielowicz [T] have 
shown how the assumption of confinement in combination with SUSY strongly constrains 
the low energy structure of the theory. The expected degrees of freedom dominating the 
low energy regime are composite operators of the gluon and gluino field which can be ar- 
ranged into a chiral superfield. These are: the gluino scalar and pseudoscalar bilinears AA, 
A75A, the corresponding gluonic quantities F 2 , FF, and the spin- 1/2 gluino-glue operator 
tr c [Fa - A]. However the program of including the purely gluonic operators ( "glueballs" ) 
as dynamical degrees of freedom turns out to be non trivial |21 01 HI IHj • In [2 El Ej the 
Veneziano- Yankielowicz low energy Lagrangian was extended so as to include all the de- 
sired low energy states, which are arranged into two Wess-Zumino supermultiplets. The 
authors of [3| pointed out on the other hand, that fulfillment of the program requires 
dynamical SUSY breaking and its consequent absence from the particle spectrum. In a 
situation where the theoretical framework seems to be still unsettled, a first-principles 
approach is welcome. This can be provided by lattice computations. 

Our goal is to verify the low energy spectrum of SYM in the case of SU(2) gauge group 
by numerical techniques. By doing this we continue past projects, see |H] for a review. 
The direct approach to the spectrum of SYM consists in studying the time-dependence of 
correlators of operators having the expected quantum numbers of the low-lying particles. 
The simplest operators of this type are the glueball, gluino-glue and mesonic operators 
also entering the low energy Lagrangians. Since gluino bilinears and glueball operators of 
the same parity carry the same (conserved) quantum numbers of the theory, it is natural 
to expect mixing among them |2j. We have to stress here that when the dynamics of 
the gluinos is taken into account beyond the valence picture, the disentanglement of the 
"unmixed" states with identical quantum numbers is not possible: only the mixed physical 
states can be the object of investigation. 1 The result is the determination of the mass 
of the lightest particle with the same quantum numbers of the projecting operator: from 
this point of view glueball and mesonic operators are equivalent. 

The action is discretized in the Wilson fashion [7] 2 where, however, the gluino is a 

1 In order to avoid confusion with the mass pattern of QCD we refrain to associate any name to the 
particle states of SYM and will refer to them according to their quantum numbers (spin and parity) . 
2 First simulations with domain wall fermions were performed in |3|. 
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Majorana spinor in the adjoint representation: 



S = S G [U] + Sf[U, A, A] ; (1) 
Sg[U] is the usual plaquette action and 

SfP, A, A] = \Y, ~ X ^ x ) ~ \ E + AW*)(r + 7m )A(x) 

X X /I 

+A(x)Vf(x)(r-7,)A(x + /i)] ; (2) 

r is the Wilson parameter set to r = 1 in our case. The gluino field satisfies the Majorana 
condition 

A = A c = CX T , (3) 

where the charge conjugation in the spinorial representation is C = 7072; the gauge link 
in the adjoint representation reads: 

[V,(x)U = 2tr[Ul(x)T a U,(x)T b ] = [V;{x)} ab = [V^x)}~J , (4) 

where T a are the generators of the color group. 

The dynamics of the gluinos is included by adopting the Two-Step Multi-Bosonic 
Algorithm (TSMB) for dynamical fermions The algorithm has the nice feature of 
accommodating any, even fractional, number of flavors. This is required for SYM since, 
schematically, the gluino has only half of the degrees of freedom of a Dirac fermion and 
consequently the fermion measure contains the square root of the fermion determinant: 
this corresponds to Nf = 1/2. In addition (cf. [TU] for details) the design of TSMB is op- 
timized to deal with light fermionic degrees of freedom, a critical factor when approaching 
the SUSY limit. Tests of the algorithm performance in QCD for light quark masses can 
be found in [TT] . 

In the Wilson discretization SUSY is broken in a two-fold way: explicitly by the Wilson 
term ensuring the correct balance between fermionic and bosonic degrees of freedom in 
the continuum limit, and softly by the gluino mass term. On the basis of the Ward 
identities H21 HB] , SUSY is expected to be recovered in the continuum limit by tuning 
the gluino mass to zero. (The situation is perfectly analogous to that of QCD, where 
chirality is recovered by tuning the quark mass to zero). However, 0(a) and 0(m§) 
SUSY violating effects are expected to distort the SUSY pattern in practical situations. 
A systematic analytical expansion in the gluino mass is missing in SYM, therefore it is not 
obvious how to set the scale for the 0{rrig) breaking (something analogous to A x = A7rf n 
in chiral perturbation theory). The only possibility, at least for the moment, to gain 
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some information on the effective "heaviness" of the gluino is to force analogy with QCD. 
Needless to say, this procedure is only of heuristic value. The strategy we adopt is to 
gradually increase the hopping parameter k in the Wilson action at fixed value of the 
gauge coupling (3 = 2.3 corresponding to a fairly small lattice size in QCD units (a ~ 0.06 
fm), pushing the simulation towards a lighter and lighter gluino. 

First large scale simulations of SYM were performed in ^0] on a 12 3 -24 for k = 0.1925. 
New sets of configurations were produced in for K = 0.194, 0.1955. We now turn to a 
16 3 • 32 lattice, whose larger time extention allows for a better analysis of the spectrum. 
We consider here k = 0.194 (simulations at k — 0.1955, 0.196 are in progress). The larger 
space extention allows us to monitor finite-volume effects in the spectrum. 

The spectrum of SYM is challenging from the point of view of numerical analysis. The 
signal for the correlators of purely gluonic operators vanishes very rapidly (ideally one 
should use anisotropic lattices). The mixed gluonic-fermionic operators, typical for SUSY 
models, receive substantial fluctuations from the gluonic content. A better asymptotic 
behavior of the effective mass, however, can be obtained by combined smearing of the 
fermionic and gluonic degrees of freedom. Finally for mesonic operators, special techniques 
are required for the disconnected term in the correlator. Here we employ stochastic 
estimators (SET) [Hj. Also, we apply an improved version [TH] of the Volume Source 
Technique (VST) [TE] . For fermions in the real representation of the gauge group, as is 
the case for SYM, the original formulation in |T6] cannot be used. The two independent 
techniques were tested in a comparative study for SYM in |15j . 

Description and some results of this study have been reported in |17j . 

The plan of the paper is as follows. In Section 2 we report details of the simulations and 
characterize the gauge sample, using analogy with QCD, by the Sommer scale parameter 
Tq and the pseudo-pion mass; the gluino mass is obtained from the soft-breaking term in 
the SUSY Ward identities; Section 3 contains methodology and results for the spectrum; 
in Section 4 we discuss results and indicate possible directions of improvement. 

2 The gauge sample 

The gauge configurations were generated by the Two-Step Multi-Bosonic Algorithm for 
dynamical gluinos jHj- We refer to JUj for a description of the algorithm. Tabled reports 
an overview of the j3 = 2.3 ensembles used in this work; the set on the 16 3 • 32 lattice 
with k = 0.194 was newly generated. The setup of the TSMB was as follows. The 
local part of the updating procedure (one cycle) consisted of two steps of heat bath for 
the bosonic fields followed by two steps of over-relaxation; the updating for the gauge 
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Table 1: Overview of the ensembles used in this work with determination of Sommer 
scale parameter and string tension. The fourth and fifth column respectively report the 
total number of configurations and of update cycles at equilibrium, the sixth column the 
number of replica lattices. 
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Figure 1: Static potential between heavy sources in the fundamental representation on 
the 16 3 • 32 lattice, k = 0.194; the line is the fit with lattice formulae. 

sector was obtained by 36 Metropolis sweeps. At the end of each cycle an accept-reject 
test was performed on the gauge configuration, along the lines of the general procedure 
described in jTOJ. Every five cycles a global heat-bath step was applied on the bosonic 
fields. The typical condition number of the squared hermitian fermion matrix was ~ 10 4 . 
The integrated autocorrelation time for the smallest eigenvalue was ~ 240 cycles. 

2.1 Static potential, string tension and Sommer scale parameter 

We measured the potential between heavy sources in the fundamental representation. 
The results on the larger 16 3 • 32 lattice confirm the picture of confinement found in [10], 
see Fig. ^ The Sommer scale parameter ro and the string tension yfa were measured by 
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Figure 2: The gluino mass as obtained from the SUSY Ward identities as a function of the 
time separation between current and insertion operator on the 16 3 • 32 lattice, k = 0.194. 
The lines indicate bounds of the fit. 



fitting the potential with the lattice formula [1 



V(t) = V + or - Arte I cos^r) 

'^(27r)3 4 V: =1 sin 2 (V2) 



r is given by 



- m/ 1 ^. w 



The results are reported in Table ^ 

Comparing the results on the two lattices, no finite-size effect beyond statistical un- 
certainty is visible in the Sommer scale parameter and the string-tension. Similarly to 
QCD, the Sommer scale parameter displays a sizeable gluino-mass dependence. A linear 
extrapolation to zero gluino mass is performed in the next subsection. 



2.2 Massless gluino limit 

In QCD the massless quark limit can be determined by inspection of the pion mass or use 
of the chiral Ward identities. In contrast, in SYM the U(l) chiral symmetry is anomalous 
and the particle with the quantum numbers of the chiral current (namely the pseudoscalar 
particle) picks up a mass by the anomaly. However, theoretical arguments pQ (cf. also 
the discussion in [T2] ) support the picture that the anomaly is originated by OZI-rule 
violating diagrams, while the remaining ones determine spontaneous breaking of the chiral 
symmetry. The diagrams of the pseudoscalar correlator respecting the OZI-rule give rise to 
the connected (one loop) term, corresponding in QCD to the pion-correlator. The analogy 
with QCD suggests the name "adjoint-pion" (a — tt) for the associated pseudoparticle: 
in the above picture this is expected to be a soft-mode of the theory, the corresponding 
mass disappearing for rrig — > 0. 
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Table 2: Quantities determined in this and previous studies: the adjoint-pion mass, gluino 
mass from SUSY Ward identities (with local SUSY current, insertion operator x , cf- 
Table 5 in [13]), spin-1/2, + and 0~ bound state masses. 

The gluino mass can be directly determined by studying the lattice SUSY Ward iden- 
tities El Il3j . where the former enters the soft breaking term. We refer to ^3] for 
the illustration of the method and discussion of theoretical aspects. One can determine 
the combination aZg m$ where Zs is the renormalization constant of the SUSY current, 
which is expected to be a (finite) function of the gauge coupling only. This quantity 
was determined in ^3] for the 12 3 • 24 lattice. We repeat here the computation for the 
16 3 -32 lattice. The results are reported in Fig. [2] where the gluino mass is plotted against 
the time separation between current and insertion operator in the SUSY Ward identities. 
Compared to the 12 3 • 24 case of 13:, the plateau establishes for larger values of the 
time separation (five compared to three); unfortunately, at these time separations the 
quality of the signal is already quite deteriorated. Table [2] contains the determinations of 
aZg rrig and am a - v in present and past works. Comparison of 12 3 ■ 24 and 16 3 ■ 32 results 
at k — 0.194 reveals a sizeable finite volume effect for the adjoint-pion mass. It should be 
noted that the sign of the effect is opposite to the usual one (however this is no physical 
mass). The gluino mass comes in larger on the larger lattice, however within a l-a effect. 

In Fig. |3] aZg nig is shown together with the squared adjoint-pion mass. The two 
quantities appear to vanish for a common value of K = k c . The estimate of k c from 
the SUSY Ward identity gluino mass, n c ~ 0.1965 [13], is not changed by the inclusion 
of the point on the larger lattice. Using this value of k c we can now extrapolate the 
Sommer scale parameter in Table^to the massless gluino situation. A linear extrapolation 
results in ro/a(m g = 0) = 8.4(4); the error takes into account the uncertainty in the 
determination of k c (assumed to be in the region k—0. 1965-0. 1975 J3])- The Sommer 
scale parameter signals the degree of "smoothness" (or "coarseness") of the gauge sample. 
In QCD, the present value would correspond to a ~ 0.06 fm (3.3 GeV), a fairly fine lattice. 
Further, assuming that the adjoint-pion drives the low energy features of SYM, as the 
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Figure 3: The gluino mass from the SUSY Ward identities and the squared adjoint-pion 
mass m a _ T as a function of 1/k (from the present and past studies [101112]); aZg on 
the 12 3 ■ 24 lattice (triangles), the same quantity on the 16 3 • 32 lattice (star); squared 
adjoint-pion mass on 8 3 ■ 16 (circle), 12 3 • 24 (boxes), and 16 3 • 32 lattice (diamond). The 
burst indicates the extrapolated massless limit from the two lightest gluino masses. 



pion does in QCD, one can estimate the degree of soft-breaking of SUSY by considering 
the dimensionless quantity M r = (rria^^ro) 2 . In QCD, validity of NLO chiral perturbation 
theory requires |21J a M r < 0.8 (corresponding to m ud < 1/4 m s ). In our case we have 
M r (n = 0.194) m 16 and M r (n = 0.1955) ~ 4.5; our lightest case would correspond in 
QCD to m u d ~ 1.5 m s . Alternatively one can consider the gluino mass from the SUSY 
Ward identity neglecting 0(1) renormalizations, again fixing the scale by the Sommer 
parameter with QCD units. In this case we obtain for our lightest gluino m- g ~ 174 MeV 
in rough agreement with the previous estimate. Since QCD and SU(2) SYM are different 
theories, the above indications are of course of qualitative nature. On the other hand, 
the relatively large average condition numbers of the fermion matrix, ~ 10 4 for k = 0.194 
and ~ 3.6 10 4 for k = 0.1955, point towards a lighter gluino. 



3 The spectrum 

As explained above, we concentrate our analysis of the spectrum on particles with spin=0 
(both parities) and spin=l/2. We investigate the glueball operators, the gluino scalar and 
pseudoscalar bilinears (meson-type operators) and the gluino-glue operator. 
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Figure 4: The disconnected component of the pseudoscalar meson correlator after sub- 
traction of the constant term on the 16 3 • 32 lattice, k = 0.194 (SET). 

3.1 Spin-1/2 bound states 

We adopt here a lattice version [T3] of the gluino-glue operator tr c [FcxA] where the field- 
strength tensor F^ v (x) is replaced by the clover-plaquette operator P^ u (x): 

0$ a (x) = ^otftr e [P ij (x)\f > (x)} ; (7) 

i<j 

only spatial indices are taken into account in order to avoid links in the time-direction. 
The clover-plaquette operator is defined to be 



with 



U$(?) = Ut{x)Ul{x + v)U v {x + fi)U,{x) 

Uj$(x) = Ul{x)U v (x - fi + ftU^x - fi)Ut(x - fi) 

u fX x ) = U u (x - v)U v {x - v - ft)Ul(x - v - p)Ul(x - p) 

U$(x) = U„(x - jl)Ut(x - fi)Ul(x + v - pL)U v {x) . (9) 

The choice of the clover plaquette vs. the regular plaquette as the gluonic field-strength 
operator in (|ZJ) is motivated by the correct behavior under parity and time reversal trans- 
formations as opposed to simply U^x). Because of the spinorial character of the gluino- 
glue, the correlator Cg a (t) has a specific structure in Dirac space. On the basis of the 
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Figure 5: Ratio of the disconnected (after subtraction) and connected components of the 
pseudoscalar meson correlator on the 16 3 ■ 32 lattice, k = 0.194 (SET). 

symmetries of the theory one can show ^2] that only two components are linearly inde- 
pendent, tr D [Cgg(x)} and tT£>[yoCgg(x)]. In our experience, the latter gives the best signal. 
In order to get a better overlap with the ground state, we apply APE smearing [22] on 
the link- variables and Jacobi smearing J22] on the gluino field simultaneously. 



3.2 bound states 

The meson-type correlators require a separate discussion because of the disconnected 
contribution. In the case of SYM one has (with A the gluino propagator): 

Cmeson(Xo ~ 2/0 ) = C conn (xo — Uq) — Cdisc{xQ ~ Vo) = 

y HTA XiV TA yiX }) (tr[TA X}X }tT[rA y J) , (10) 

X X 

with r G {1, 70} (observe the factor 1/2 reflecting the Majorana nature of the gluino). The 
disconnected term requires the estimation of the time-slice sum of the gluino propagator 

S a p{xo) = 5Jtr c [A xcw g] . (11) 

x 

For this, we use the stochastic estimator technique (SET) [2] with complex Z 2 noise 
in the spin explicit variant SEM ^H]- In this case each estimate of the time-slice sum 
is obtained by inverting the fermion-matrix with source (wg ) X bp — $apVxb wnere Vx°b 
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the last case only a rough indication of the first excited state mass could be obtained. 
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are independent stochastic variables chosen at random from ^=(±1 ± i). Here we use 
point-like operators (i.e. no smearing on the gluino). 

On the larger 16 3 • 32 lattice the computed meson correlator displays an offset: its 
long-time behavior is not purely exponential, since a constant term also appears. Such 
a constant term is theoretically excluded in the correlator by the symmetries of the the- 
ory. It is present in both SET and VST (see below) determinations of the disconnected 
contribution and does not decrease by increasing the number of the random estimators. 
In contrast, it is absent in the connected correlator. We conclude that its origin is to be 
traced to some cumulative numerical effect in the stochastic computation of the discon- 
nected contribution. 

In Fig. |3] we show the disconnected component of the pseudoscalar meson correlator 
after subtraction of the constant term. In Fig. El the ratio between the subtracted discon- 
nected component and the connected one is reported. For a comparison with the same 
quantity in the case of QCD see e.g. |2HES1; we remark here that the case of SYM is quite 
different since the connected correlator is not related to a physical particle, but rather to 
the pseudoparticle a — n discussed above. 

We cross-check the SET with the improved version ^3] of the Volume Source Tech- 
nique (VST) ^B|, applying to fermions in real representations of the color group. The 
improvement consists in averaging the time-slice sums over random gauge transforma- 
tions and therefore eliminating the gauge non-invariant spurious terms. The two methods 
deliver consistent results of comparable quality at similar computational cost. For the 
sake of brevity, we present here only those from SET. 

Another operator with the right quantum numbers (0~) is the pseudoscalar glueball 
operator. This is given by a linear combination of closed loops of link variables which 
cannot be rotated into their mirror image (cf. e.g. |10|). We considered the simplest loops 
of this kind. Unfortunately this operator does not give a clear signal on our samples. 

3.3 + bound states 

Since the meson-type correlator does not show any appreciable signal for the + state, we 
turn to the scalar glueball operator. The standard operator in this case is 

O glueball (x) = tr c [U 12 (x) + U 23 (x) + U 31 (x)} . (12) 

We use fuzzy operators by applying APE smearing on the link variables. 
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Figure 9: Mass of lightest bound states of SYM determined in this work. The shaded 
region represents the presumed location of the massless gluino on the basis of the SUSY 
Ward identity analysis. 

3.4 Results 

For all particles we measured the effective masses (Figs. HUB))- in many cases a clear 
plateau could not be determined. In order to get a better determination of the ground 
state mass we used constrained two-mass fits (bounds in the figures). In some cases (the 
spin-1/2 particle on the larger lattice at n = 0.194, Fig. and the pseudoscalar particle 
at the same n value, Fig. |7|) the two-mass fit can be cross-checked against a plateau of the 
effective mass. We ensured the stability of the two-mass fits by systematically varying fit 
ranges (for details see The effective mass of the pseudoscalar meson on the 16 3 ■ 32 

lattice, lower panel of Fig. [7[ was determined after subtraction of the constant term in 
the correlator discussed in Subsection 13.21 

In the case of the scalar glueball operator, a decrease of the signal/noise ratio was 
observed on the larger lattice, as a consequence of which no determination of the mass 
was possible. 

Results on the determinations of the ground state masses are reported in Table 121 and 
Fig El A discussion of the results will be presented in the following section. 



4 Discussion 

Our analysis of the low-lying spectrum of SYM shows a slow approach of the correlators 
to the asymptotic behavior where only the ground state dominates. This is evident in the 
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case of the gluino-glue and mesonic correlators; in the case of the glueball correlator, the 
quality of the signal is not good enough to make definite statements. Excited states with 
masses comparable to that of the ground state are strongly coupled to these operators. In 
some cases a plateau of the effective mass emerged and consequently allowed us to cross- 
check the results of the two-mass fits. The situation can be improved by implementing 
optimized smearing on the operators (in the case of the gluino bilinears we use point-like 
operators) . 

The excited states which hamper the determination of the ground states are of physical 
interest by themselves. According to j2j EJ |H] , the first excited states should be arranged 
in a second Wess-Zumino supermultiplet. The "higher masses" in our two-mass fits can 
give a first indication of the masses of these excited states: the ground state masses lie 
in the region 0.2-0.5 (in lattice units), while the higher masses are in the region 0.8-1. 
A more refined analysis of the excited states, however, could be obtained with matrix 
correlators. In the scalar sector, one would naturally include the gluino scalar bilinear in 
addition to the glueball operator. Given the large fluctuations observed on the former, 
the employment of variational methods would then be advisable. A similar analysis 
could be done in the pseudoscalar sector with the corresponding gluino bilinear and the 
pseudoscalar glueball operator. 

A more fundamental question is whether the employed operators are optimal in the 
sense of maximal overlap with the low-lying bound states of SYM. Investigations could 
go in the direction of different operators and different quantum numbers [2*5] . 

In the following we restrict the discussion of our results to the ground states (Fig. [""J). 
One of the goals of this study was to check finite volume effects by comparing lattices with 
different spatial extension. This can be done for our value of k = 0.194 where two different 
lattice sizes are available, L s = 12 and 16. The direct comparison shows, see Table [*"J that 
a sizeable deviation is present for the pseudoscalar particle. Contrary to expectations, 
the particle comes in heavier on the larger lattice. For this lattice however an unexpected 
constant term is observed in the long-time behavior of the correlator, which could hint 
at some systematic effect in the stochastic determination of the disconnected correlator 
on large lattices. The pseudoscalar particle is the lightest particle for our lightest gluino 
(« = 0.1955), though in this case only data for the smaller lattice is available. The scalar 
and the spin- 1/2 particle have comparable masses, compatible within errors. 

Conclusions on the relevance of soft breaking terms require the control of finite lattice- 
spacing effects. Using analogy with QCD in absence of other indications, we argue that 
our mesh is relatively fine, while the gluino is still quite heavy. Next steps will be therefore 
to consider larger values of k on large lattices. 
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